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ABSTRACT 


The Naval Postgraduate School (NPS) Micro-Air Vehicle (MAV) wind tunnel 
was refurbished in this study. This wind tunnel has a 61 by 38 centimeter test section 
with a contraction ratio of 6.75, and testing speeds up to 9 meters per second (m/s). The 
objectives of this work were to create a high quality, customized facility for further MAV 
study, and to fully characterize the wake of the MAV at a plane downstream of the 
model. Extensive repairs were made to the wind tunnel intake, test section, and fan. The 
flow field produced in the tunnel was analyzed with hot wire anemometry. The 
turbulence intensity ranged from .7% to .5% for freestream velocities of 2 m/s to 5 m/s. 
The velocity variation across the test section ranged from .251% to .125%. A model of 
the NPS MAV was placed in the test section and the wake was studied with and without 
the main wing, and in vertical and horizontal orientations. The wake was characterized 
with hot wire anemometry and flow visualization techniques at various flapping 


frequencies and freestream velocities. 
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I. INTRODUCTION 


A. BACKGROUND 


In the past, most military airborne surveillance missions were performed by 
expensive large piloted aircraft. Recent technological advances have led to the use of 
Unmanned Aerial Vehicles (UAVs) for reconnaissance tasks. UAVs are advantageous 
because they remove the human from the battlefield, are much smaller and less 
expensive, and are much stealthier than typical piloted aircraft. The next level of 


development will be in the area of Micro Air Vehicles (MAVs). 


Real time satellite data, spy plane and UAV surveillance are becoming more 
accessible. However, the information provided from these sources is frequently not 
tailored to the small area detail required by lower level military leadership. The squad 
leader must still put troops in harm’s way to scout over the next hill, through the alley or 
inside the next building. An effort is underway to change this. In 2000, the Defense 
Advanced Research Projects Agency (DARPA) set guidelines for the first generation of 
MAVs. Since then, major MAV competitions have taken place on national and 
international levels. Specifications and mission requirements vary, but the trend is clear: 


MAVs will play a significant role in future military action. 


MAVs are especially advantageous in missions where human risk is involved. 
They could be used by government agencies involved in terrorist or hostage situations, 
search and rescue operations, environmental monitoring of smokestack emissions or 
chemical spills. Military missions such as video or infrared surveillance, battle damage 
assessment, weapon targeting, communications relay, detection of toxic gasses, radiation, 


or detection of land mines are also suited for MAV use. 


B. OBJECTIVE 


This project chronicles the rebuilding efforts of the Naval Postgraduate School 


(NPS) Micro Air Vehicle (MAV) wind tunnel, which will be used for MAV testing and 
1 


research. Extensive repairs were made resulting in a facility far superior to that used 
previously. The NPS low speed wind tunnel formerly used for MAV testing has been 
decommissioned and slated for removal. This facility had a number of flaws including a 
poor lip design and a close proximity to building walls and other articles that led to a 
turbulence intensity in excess of 2% at MAV velocities. It suffered from the ill effects of 
a short, rapid 9:1 ratio contraction and only one layer of screen in the intake with 
insufficient settling area. The details of the detriments of this construction will be 


developed further in subsequent chapters. 


The goal of this work was to provide a viable facility for MAV testing and 
research. Additionally, this paper sought to analyze the flow around the NPS MAV using 


Constant Temperature Anemometry and flow visualization. 


C; MAV CHALLENGES 


iL, Low Reynolds Number Flight 

Because of their small size and low airspeed, MAVs operate at Reynolds numbers 
(Re) lower than conventional craft. One cannot simply scale down an existing aircraft to 
a six inch size and expect it to fly. Aerodynamic behavior is different in this low Re 
regime and has not been studied in depth. Most aircraft fly at Re in the millions, even 
100 million for the largest commercial jets. At Re above 200,000 many simplifications 
can be made to the governing equations of flight to facilitate the design process. These 
simplifications are allowable because the flow around airfoils in this regime is dominated 
by inertial forces and can be modeled with vortex panel codes. At low Re, viscous 
effects play a larger role and can not be ignored. This required the use of Navier-Stokes 
based codes, which are more accurate, but also more time consuming and expensive. The 
NPS MAV operates at around Re 20,000 for the main wing and 5,000 for the flapping 


wings requiring these Navier-Stokes based methods. 


Airfoil performance in the low Re regime is sensitive to separation. If the flow 
becomes separated from the airfoil, it will lose lift in a condition known as stall. Airfoils 


develop a pressure profile based on their shape and angle of attack. In general, the lowest 
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pressure develops on the top of an airfoil at a position 30% to 50% of the chord length 
known as the center of pressure (Cp). From the leading edge of the airfoil to this point 
the pressure gradient is from high pressure to low pressure known as favorable, or acting 
in the natural direction of the flow as the wing moves through the air. From the center of 
pressure to the trailing edge the gradient is unfavorable, or adverse, leading to potentially 
serious problems. An adverse pressure gradient tends to slow the flow as it moves from 
low pressure to high. Near the surface of the wing, the boundary layer has already 
reduced the momentum of the flow, so the adverse pressure gradient exacerbates the 
issue. The flow will gradually slow to zero velocity and even reverse direction. This 
recirculation flow is usually not attached to the surface of the wing and the flow is 
separated, or the airfoil is stalled. When stalled, the airfoil has dramatically reduced lift 
and the craft quickly loses altitude. This is a particularly dangerous condition for low 
flying vehicles like MAVs. At higher angles of attack, the pressure at the Cp decreases. 
This increases the adverse pressure gradient to the trailing edge and makes stall more 
likely. The NPS MAV handles the problem of stall with its unique design described in 
Part D below. 

Zi Additional Challenges 

There are many other engineering challenges concerning possible MAV 
requirements. The vehicle must be lightweight and compact with flight stabilization 
necessary for low Re flow as outlined above. This requires the use of high energy 
density propulsion systems and high strength to weight ratio structural components. 
Lightweight, low power on board communication equipment must be integrated into the 
design with sufficient bandwidth for real time imaging. On board guidance, navigation 
and geolocation systems may also be desirable. The development of the sensors 


themselves will likely prove difficult as well. 


D. THE NPS MAV 


To meet these challenges the NPS MAV has taken an unconventional biomorphic 
approach. The craft consists of low aspect ratio wings flanking the fuselage followed by 


a pair of identical counterphase flapping wings. This allows the flapping wings to take 
3 


advantage of a phenomenon known as ground effect for increased performance. Ground 
effect is caused by the interruption of wingtip vortices. These vortices destroy large 
amounts of lift generated by the wing through increased downwash behind it. The 
symmetric action of the flapping wings interrupts these lift destroying vortices [REF 13]. 
This dual flapping wing design is also mechanically and aerodynamically balanced. A 
single wing system would be unbalanced and less efficient as the net torque on the main 
structure would induce undesired motion on the fuselage and main wing. Looking at the 
flow 2 dimensionally, a flapping wing sheds vortices at its tip adversely affecting its 
performance. However, for a pair of flapping wings, the image wing is positively 


affected by phenomenon enhancing the overall systems performance. 


Another aerodynamic phenomenon the NPS MAV is uniquely suited to handle is 
the flow separation and stall problem at low Re flight mentioned above. The flapping 
wings produce a favorable pressure gradient in front of them and on the adjacent aft 
portion of the main wing. This suppresses flow separation and stall even in the low Re 


flight regime used by MAVs [REF 13]. 


The NPS MAV weighs around 12 grams and currently has no payload. It can fly 
between 2 and 5 m/s for a great degree of maneuverability. It currently does not have 


active pitch control, limiting the maximum safe speed. 


Il. RENOVATION 


A. OVERALL DESCRIPTION OF THE FLAWS 


The wind tunnel that would become the Naval Postgraduate School Micro Air 


Vehicle tunnel was discovered in a basement storage area as shown in figures | and 2. 





Figure 1. | Wind Tunnel in Storage 





Figure 2. Exhaust Section in Storage 


It was mostly complete but damaged in several sections and covered in a thick 
layer of dust. The dust itself was a flaw as dust always accumulates on the honeycomb 
and screens in irregular patterns which alters the flow [REF 1]. The honeycomb flow 
straightener material was damaged in several places, including a serious defect in the 
lower left corner of the intake (see Figure 3). The honeycomb was constructed of two 
separate sections. Several of the cells were blocked along the joint right through the 
center of the tunnel by the epoxy used to attach the two parts (See Figure 4). These 


blocked cells would have led to a disturbance in the test section flow [REF 1]. 
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There were numerous holes in various parts of the test section including a large 
hole in the side of the test section that was inaccurately sealed with a plug. There was a 
felt like material glued to the bottom of the test section. These defects are shown in 


Figure 4. All of these defects would have certainly disturbed the flow in the test section. 





Figure 5. | Test Section Access Hole and Felt on Floor Section 


The motor used to drive the fan was a 1.5 hp General Electric AC motor. This 
device had limited speed control severely restricting its utility. To alter the tunnel flow 
speed the operator had to secure the flow, remove the stator, loosen the outer fan hub and 
adjust the fan blades before reassembling and checking the flow. The blades would have 
to be carefully aligned to the same pitch for aerodynamic balance and an even test section 
flow. This was an inaccurate and tedious method of adjusting the flow that proved 
unacceptable. More importantly, in short wind tunnels with relatively small contraction 
ratios such as this, the fan speed should be varied rather than the blade pitch to maintain 


an even velocity front in the test section [REF 1]. 


The tunnel also experienced significant vibration during operation that would 
limit its usefulness in MAV testing. This was caused by an unbalanced fan shaft. 
Vibrations must be limited because they adversely affect measuring systems, model 
operation, and flow visualization. Furthermore, vibrations can lead to transition and 


separated boundary layers. 


A suitable location for the wind tunnel was found. When attempting to move the 
tunnel, it was discovered that the wheels had sunk into the asphalt floor. The intake 
section was raised with a hydraulic jack and plywood inserted under the wheels. 
Measurements were taken and it was discovered that the tunnel would not fit along the 
path to the new location. To accommodate the move, the intake section had to be 


disassembled. With this complete, the facility was moved to the new location. 


B. WIND TUNNEL INTAKE 


The wind tunnel intake was completely disassembled to facilitate moving the 
structure. The intake lip was attached to the structure with thirty-six bolts. The 
honeycomb and the layers of screen were held in place in individual wooden frames. The 
three frames were held together by thirty-two threaded rods with nuts on each end. The 
honeycomb flow straightener was judged to be beyond repair. It was carefully cut from 
its frame. The remaining epoxy was removed and the frame prepared with a belt sander. 
A replacement material was obtained but delivered with some damage. The corners were 
collapsed and a few areas in the center were slightly deformed. A heavy gauge wire was 
inserted into the cells to carefully reshape the material. The new honeycomb material 
was slightly collapsed with fine aluminum sheets to fit inside the frame. Then the 
honeycomb was permanently fixed to the frame by painstakingly applying epoxy around 


the interface. 


One layer of the original screen was wired together from two separate sections of 
material and needed replacement as shown in Figure 6. This represented an improper 
splice that would have introduced turbulence into the flow [REF 1]. When an intake 


section is too large for a single section of screen, splicing must be done very carefully by 
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brazing individual wires together. When suitable replacement materials were obtained; 
the screen was replaced by pulling tight and stapling a single piece of appropriately sized 


screen over the frame. 





Figure 6. Damaged Screen Section 


C. TEST SECTION 


The test section ceiling was mounted on failing hinges. New ones were obtained, 
holes drilled and tapped, and the top section was firmly attached to the rest of the 


structure as shown in Figure 7. 
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Figure 7. | New Hinges 


The felt material on the test section floor was removed with a commercial solvent 
and the area is shown clean in Figure 8. The large hole in the test section also required 
repair. New sheets of optical grade acrylic were obtained. This material was chosen to 
allow Laser Doppler Velocimetry measurements in future application. A_ large 
rectangular shape was removed from around the hole in the test section with a jigsaw as 


shown in Figure 8. 





2 = 
Figure 8. | Test Section Cut for Optical Quality Upgrade 
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The area was carefully sanded until the new sheets could be inserted with a tight 
fit. Small strips of aluminum were obtained to support the pieces. Holes were drilled and 
tapped for mounting the aluminum strip to the tunnel and to the new sheet of material. It 
was cut to final size and a slot was added for flow measurement. Additional sheets have 
been prepared for easy replacement if necessary. A thin weatherstripping material was 
added at the joints to limit leakage. Large holes were cut in the framing of the tunnel to 


accommodate lighting or equipment mounting from the floor as shown in Figure 9. 





Figure 9. _—_ Holes in Frame for Floor Mounting 


D. FAN 


The tunnel vibrated excessively when operated. It was determined that the shaft 
and hub were out of balance. The fan section was removed for disassembly. Simple 
static balancing showed the structure to be significantly unbalanced. There were several 
holes drilled in the hub which contributed to the unbalanced condition as they appeared to 
be drilled on the light side of the hub (Figure 10). Three holes were drilled on the heavy 


side but were not having a significant effect. Weights were attached to the blades to 
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achieve balance and it was determined that simply drilling holes would not be enough to 
correct the situation. Some of the existing holes were filled with specially fabricated 
steel inserts and a weighted washer was constructed to bring the assembly in balance 


(Figure 11). 





Figure 10. Unbalanced Fan Hub 
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Figure 11. Balanced Shaft 


The fan was remounted with the tunnel completely assembled and tested. The fan 
stator would no longer fit with the addition of the balancing weight. One quarter inch of 
material was removed from the stator with a rotary tool and a cutoff wheel to allow full 
reassembly. Upon testing, the vibration was no longer present. However, while running 
at high speed the fan exhibited a grinding noise and the test run was halted. The shaft 
was difficult to turn by hand and the fan section was once again removed. It was 
determined that the bearings were failing. The shaft was supported by two ball bearing 
units under the hub and a flange bearing at the bottom of the fan section. In order to split 
the fan section the lower flange bearing had to be detached. The bearing was frozen on 
the shaft and proved difficult to remove. A specially modified two jaw gear puller was 
obtained. The ends of the extracting arms were sharpened on a grinder to fit in the small 
space between the bearing and the lower fan section. With extreme care, the tool 
successfully removed the bearing. This allowed separation of the lower fan section, after 
removing the bolts and prying it apart. At this point, the fan hub and blades required 
removal. The blades and corresponding sections of both halves of the hub were 
numbered to ensure the same balance upon installation, and the blade angles were marked 


to ensure the same aerodynamic balance. The bolts were removed and the modified gear 
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puller placed back in service to remove the outer fan hub section. With this removed, the 
blades were freely extracted. An additional gear puller was used to remove the inner fan 
hub and locating Woodruff key. The shaft was then extracted from the bottom of the 
upper fan section leaving the structure shown in Figure 12. This exposed the two ball 
bearing units, one on each end of the support shown in the middle of the photo. They 


proved difficult to remove as well. 





Figure 12. Upper Fan Section 


To remove and install the new ball bearing units, a pulling system was devised. 
This consisted of thick scrap aluminum, mill mounting pieces and a threaded rod. A 
piece of rectangular scrap was inserted into the bearing hole and righted. The threaded 
rod was then screwed into this piece which would apply force to the inside of the bearing. 
A thick section of scrap with a through hole for the threaded rod was inserted and 
balanced against the structure as shown in Figure 13. A nut was then threaded on the rod 
and turned to apply force to the inside of the bearing for removal. The bearing unit on 


the other end of the support was removed in this manner as well. 
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Figure 13. Bearing Removal System 


Installation required similar tools. After new bearings were obtained, a longer rod 
was inserted through the whole support structure and threaded into a scrap piece on the 
other end. Then the same setup shown in Figure 13 was utilized except the support 
pieces contacted the new bearing. The nut was tightened and the bearing was pressed 


into place. The same action on the other end installed the opposing bearing. 


The fan shaft was found to be corroding in several places. For this reason, it was 
placed in a lathe and polished for easier reassembly. The upper fan housing mating 
surface was left with epoxy residue as shown in Figure 12 above. This surface and the 
mating surface on the lower fan housing were prepared with a belt sander. Thin strips of 
weather-stripping material were added to ensure an air tight fit. The shaft was then 
inserted into two new bearings in the support housing. The inner fan hub was pressed 
into place and located with the key. The blades were reinstalled in their correct mounting 
holes and the outer hub was pressed into place in its corresponding position. The lower 
fan housing was then reattached and bolted together. The new flange bearing was added 


to the shaft and bolted to the lower fan housing. The coupling was added to the end of 
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the shaft and the completed fan housing was reattached to the wind tunnel structure, with 
additional weather-stripping inserted in the joint. The completed assembly is shown in 


Figure 14. 





Figure 14. Fully Assembled Fan Section 


E. MOTOR 


An alternative driving force for the wind tunnel was sought to improve the control 
of the airflow. It is more desirable to change airflow velocity by changing blade speed 
than by changing blade pitch [REF 1]. A direct current motor and controller was 
recycled from a neighboring wind tunnel to serve this purpose. Unfortunately, the shaft 
did not align with the fan when the motor was mounted and it required a smaller coupling 
to accommodate the smaller shaft diameter. The mounting surface consisted of an 
aluminum plate bolted to the wooden frame of the wind tunnel. The plate was removed 
and the wooden surface was cut away with a router to allow the mounting plate to be 


recessed as shown in Figure 15. The motor was then installed and connected to the fan 
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shaft with a new flexible coupling. The system was started and it was discovered that the 
flow was in the wrong direction. Research determined that this motor could be operated 
in either direction so the controller was disassembled and rewired to allow the motor to 


provide the correct flow. 





Figure 15. Recessed Motor Mounting Plate 


F. TRAVERSE MECHANISM 


The traverse mechanism used in this project was recycled from a neighboring 
wind tunnel. It was a three axis unit of aluminum construction fabricated by Velmex and 
shown in Figure 16. In what would become the x-axis in the orientation used in this 
project, the mechanism was out of alignment. Upon further inspection, it was discovered 
that the belt driving the twin lead screws was loose. Additionally, the drive nut on one of 
the slides was loose and slipping. The mechanism was disassembled, cleaned, and 
lubricated. All of the drive nuts were tightened enough to provide a positive grip on the 
lead screw without stalling the motors. The belt shown in Figure 17 was tightened and 


the system was reassembled. 


18 





Figure 16. Traverse Mechanism 





Figure 17. Traverse Mechanism Drive Belt 


Next, a suitable mounting table was obtained. The table surface was raised to an 
appropriate height for the wind tunnel test section. Aluminum L-brackets were fabricated 
for locating the mechanism on the table. The table was then bolted to the wind tunnel 
frame for a positive and stationary location. The mechanism was shimmed to a level 
position. A bracket was then fabricated for attachment of the anemometer probe holder 


to the z-axis slider of the traverse mechanism. Then a probe holder was constructed of 
19 


carbon fiber reinforced spruce wood and was clamped to the mechanism bracket with a 
pair of C-clamps. The end of the probe holder was carefully dilled and cut to receive the 
anemometer probe and thermocouple. A thin plastic screw was utilized to clamp the 
probes in position at the end of the holder. The holder was shaped in an elliptical cross 


section to minimize wind disturbance in the flow field. 
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HI. EXPERIMENTAL SETUP 


A. WIND TUNNEL 


1. Description 

The NPS MAV wind tunnel is an open circuit design with the intake and test 
section mounted horizontally. After the test section, the flow is turned upward through 
guide vanes and driven by an axial fan powered from an electric motor. The facility has 
an overall length of 13 feet. It is 84 inches tall at the tip of the stator cone at the exhaust 
and at the upper intake lip. The entire structure is supported by a wooden frame with 
casters for mobility. The intake begins with a half section of 10 inch diameter piping 
along all four sides as a bellmouth, as is used in the University of New Mexico facility 
[REF 7]. A bellmouth type intake helps reduce losses, but, more importantly in this 
application, it helps provide a more uniform flow [REF 9]. The settling chamber consists 
of the honeycomb section along with two layers of screen and their associated settling 
regions. Honeycomb followed by screen combinations have proven to be the best 
method of reducing turbulence [REF 12]. This completes the intake and marks the 
beginning of the contraction. A splitter plate is not used in this facility. The intake and 


contraction section are illustrated in Figure 18. 


The honeycomb section is “a guiding device through which the individual air 
filaments are rendered parallel” according to Prandtl [REF 11]. The section used in this 
facility had a cell size of .25 inches and a thickness of 1.375 inches. This yields a length 
to size ratio of 5.5. Almost complete elimination of the lateral components of turbulence 
occurs in the range of 5-10 [REF 8]. The section was aluminum in construction. This is 
superior to impregnated paper honeycombs sometimes used because the aluminum 
provides more precise dimensions and leads to better flow characteristics [REF 8]. The 
cell size should be smaller than the smallest lateral wavelength of velocity variation. 
This means roughly 150 cells per settling chamber ‘diameter’ or 25,000 cells in total for 
adequate results [REF 8]. The facility has approximately 38,800 cells and the additional 


number of cells compensates for the slightly smaller than typical length to cell size ratio. 


21 


Honeycombs themselves shed small scale turbulence, so the honeycomb should be placed 


upstream of the screens to minimize the turbulence level, as is done in this tunnel. 


Screens in the settling chamber make flow velocity more uniform by imposing a 
pressure drop proportional to the velocity squared. They also tend to reduce the 
boundary layer thickness and refract incident flow toward the normal. The most 
important function in this application is the ability of screens to reduce turbulence 
intensity. The screens were made of .01 inch diameter wire with a mesh size of .05 
inches. The screens should be placed 50-75 mesh diameters apart so that the screen shed 
turbulence has decayed [REF 12]. The screens in the facility are spaced at 3.25 inches in 


their mounting frames for a spacing of 65 mesh diameters. 


The intake is 45 inches high and 54 inches wide yielding a contraction ratio of 
6.75. The main effects of a contraction section are to reduce the mean and fluctuating 
velocity variations to a smaller fraction of the average and to increase the average mean 
velocity. Contraction ratios between 6 and 10 have been found to be adequate in small, 
low speed tunnels [REF 8]. The wall shape of a contraction is based on producing a 
uniform steady stream at the outlet of the section. This would dictate a long contraction. 
However, this must be balanced with the desire to minimize the size of the boundary 
layer, favoring a short contraction. A short contraction section also moves the flow 
closer to boundary layer separation and flow unsteadiness. Additionally, the concave 
shape at the contraction inlet has a destabilizing effect on the boundary layer while the 
convex shape at the exit tends to stabilize the flow. A desirable design would then dictate 
that separation is just avoided (a minimum length) and the nonuniformity at the exit is 
tolerable, typically less than 1% variation in mean streamwise velocity outside the 
boundary layers [REF 8]. From the flow analysis on the tunnel it will be shown that the 


contraction design is satisfactory. The contraction section is shown in Figure 18. 
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Figure 18. Intake and Contraction Sections 


The length of the contracting section is 40 inches, thus the length to height ratio is 
.89, within the range of similar wind tunnels and exactly the same as the Mixing Layer 
Tunnel at NASA Ames [REF 9]. The contraction is three dimensional. The wall contour 
shape is unknown; however the statistics are similar to the fifth order polynomial 
contraction used in the boundary layer tunnel at NASA Ames [REF 9]. While the actual 
form of the contour is not that important, the smoothness in the contour shape is critical. 
Figure 18 shows that this section is quite smooth. It is also important that the ends be 
smooth and the second radius of curvature be less than the first, and that the ends join the 
parallel sections so smoothly that the first and second derivatives of the curvature are 


close to zero [REF 8] as shown in Figure 19. 
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Figure 19. Contraction Contour [From REF 14] 


The test section is 53 inches long and made of plexiglass. The test section is 
rectangular and measures 15 inches high and 24 inches wide. The cross section design is 
not critical and depends on the needs of the user [REF 8]. The entire ceiling of this 
section can be raised for access necessary for testing. There are two supports over the 


test section length for mounting to the frame. 


The airflow is driven by an axial fan near the exhaust. While prerotation vanes 
are normally used, it is acceptable for a lightly loaded fan, as in this application, to use a 
set of straightener vanes downstream of the fan as long as the number of fan blades is 
different from the number of straightener vanes [REF 9]. The fan has 4 blades while the 
stator has 9 blades and a center cone. The fan is turned by a Minarik Electric direct 


current motor and a Motor Master controller. 


B. TRAVERSE MECHANISM 


The traverse mechanism was a Velmex Unislide model B4026P40J with a Bodine 
Electric DC stepping motor type 34T1BEHD attached for motion in the x-direction. The 
y-direction was controlled by a Velmex Unislide model B4027P40J with a Superior 
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Electric stepping motor M091-FD09, and the z-axis was regulated by a Velmex Unislide 
model B2515 P40J with A Bodine Electric DC stepping motor type 23T2BEHD attached. 
All of the motors were capable of 200 steps per revolution. This allowed for 16,000 
counts per inch of movement when combined with the lead screw pitch. The mechanism 


was controlled by a Velmex VP 9003 controller. 


C. CONSTANT TEMPERATURE ANEMOMETER 


1. Equipment Used 

The constant temperature anemometer (hotwire system) used was a Dantec 
Streamline 90N10 frame with five 90C10 CTA modules. The data was processed by a 
Pentium computer running Streamware software version 3.4 after being fed through a 
National Instruments A/D board model BNC-2090. The hotwire system was calibrated 
with a dedicated calibration unit, the Dantec Streamline 90HO2 flow unit code 
9090H0021 serial 149. 

Z Theory of Operation 

A hotwire flow measurement system consists of the components shown in Figure 
20 below. The probe used was a miniature wire unit as this type is recommended for 
unidirectional gas flows with turbulence intensity up to 10% because they have the 


highest frequency response. The probe consists of a 5 44m diameter tungsten wire 1 mm 


long supported by stainless steel prongs at each end. The prongs are then coated in a 
ceramic with the ends connecting to a probe support and connected to the anemometer 


through a probe cable. 
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Figure 20. Hotwire Measuring System [from REF 10] 


The probe wire’s resistance changes with its temperature. This then becomes the 
measurement leg of a Wheatstone bridge opposite a variable resistor, which defines the 
operating resistance, and hence the operating temperature of the instrument. The system 
also contains a thermocouple in the flow field next to the hot wire probe for temperature 
compensation. When the bridge is in balance, no voltage difference exists across its 
diagonal. If the flow velocity were to increase, more heat would be transferred away 
from the wire and its temperature would go down so its resistance would go down 
resulting in an error voltage at the input of the current regulating amplifier. This will 
cause the probe current to increase. The wire will heat and increase in resistance until the 
balance is restored. The current regulating amplifier has a very high gain so a condition 
of bridge balance exists almost independent of the flow velocity past the wire. The 
response time of the wire is thus reduced to a few micro seconds. The probe current is 
calculated by the voltage drop across the bridge. In this way the system calculates flow 


velocity at very high frequencies based on the heat transfer from the wire. 


D. MAV MODEL 


The NPS MAV model was constructed by Professor Kevin Jones. It was similar 
to the flying MAV but without the rudder. The model had a much larger fuselage to 


house the more robust motor and transmission, while the flying craft places motor, 
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transmission, battery, and radio gear in a smaller space. This allowed the model to 
operate indefinitely, while the flying MAV motor has a limited lifetime. Additionally, 
the main lifting wings in the model shown below are covered with Japanese tissue paper 
while the operational aircraft used the same microfilm plastic sheeting that covers the 


flapping wings. 














Figure 21. NPS MAV Model [From REF 4] 


The frame of the main wing was composed entirely of light weight balsa wood 
and carbon fiber assembled with a cyanoacrylate adhesive in a sandwich fashion. The 
main wings were slightly wider than the flapping wings at 10.5 inches with a mean 
geometric chord of 5.5 inches. The main wings were rectangular in shape and mounted 


with 8.5 degrees of dihedral. 


The flapping wings were powered by a small electric motor run from an external 
power supply for this project. To power the model, an Elenco model MX-9300 power 
supply was used. The model flapping frequency was measured using a strobe light. A 
chart was developed relating motor voltage to flapping frequency and used during the 


experiment. 


P| 


The flapping wings were driven from the motor, through the transmission and to a 
crankshaft. The crankshaft converted the rotation of the motor to linear motion required 
for the flapping wings by driving connecting rods to each flapping beam. The flapping 
beams were of carbon fiber sandwich construction common to the other structural 
members. The flapping beams were attached to the front of the fuselage on one end and 
the flexible wing mounts on the other. The plunging motion of the flapping wings was 
caused by the motor’s rotation while the pitching motion was aeroelastic. As the flexible 
wing mounts at the end of the flapping beams flexed, the pitch was induced in the 
flapping wings. The amount of pitch was related to the stiffness of the flexible wing 
mounts. The flapping wings functioned in counterphase to eliminate the inertial and 
aerodynamic effects of a single wing moving up and down. The flapping wings were 
identical in construction, each having a chord of 1.5 inches in the midspan and tapered at 
the ends with a span of 9.75 inches. The flapping wings are tapered to limit spanwise 
oscillations of the wing and provide no noticeable decrease in thrust. The leading edge 
was constructed of balsawood dowel. The surface was covered with thin plastic film 
supported at seven points with chordwise ribs of carbon fiber sandwich construction. 
Two layers of the rib material run two thirds of the chord length and a third layer runs to 
the trailing edge of the wing. Although multiple crankshafts and connecting rods were 
available to adjust the plunging amplitude and mean separation, for the experiment, the 
plunging amplitude, from the centerline to top-dead-center, is 1 inch or 0.67 chord 


lengths, with a mean separation of 3 inches. 


E. PITOT TUBE FLOW REFERENCE 


A small pitot tube was used as a reference free stream velocity for different data 
runs. The tube was connected to a MKS Baratron type 223BD pressure transducer 
powered by an Analog Devices type 904 power supply. The voltage is read by a Metex 
ME-11 voltmeter. This system was calibrated by the hotwire and the resulting chart is 


located in Appendix B. 
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IV. EXPERIMENTAL PROCEDURE AND RESULTS 


A. CALIBRATION 


The hotwire system was calibrated using its own calibration module for maximum 
accuracy. The calibration module operates from a pressurized air supply and creates a 
free jet. The free jet is produced through a high resolution control valve and a series of 
Laval nozzles that provide a stable mass flow regardless of back pressure. Thus the free 
jet velocity is highly accurate and adjustable to fractions of a percent. The probe and 
thermocouple were attached to the module with the hotwire sensor positioned directly in 
the flow. The module was set to calibrate from | to 10 m/s in accordance with the testing 
velocities of interest in 10 points arranged logarithmically as recommended by the 


manufacturer. This resulted in the table in Appendix C with a fourth order polynomial fit 


+C E> _+C,E.,,+C,E- 


corr corr corr * 


according to the equationU =C,+C,E 


corr 


B. FLOW VERIFICATION OF THE TUNNEL 


Le The Grid 

The grid used in the flow verification of the tunnel consisted of a plane in the x or 
streamwise, direction. The y direction was taken horizontally with the z direction in the 
vertical. The origin was taken as the center of the tunnel with the negative y direction 
running to the left side of the center and the negative z direction extending below the 
center of the test section. Data was taken at 31 positions in the y direction at 7 different z 
positions. Boundary layer effects were captured on the left side starting .125 inches from 
the wall. Ten points were taken within an inch of the wall and the rest of the points were 
spaced one inch apart across the tunnel. The traverse mechanism did not have sufficient 
travel to obtain boundary layer information from both sides of the test section. It was 
assumed to be similar on all sides. No vehicle testing will ever be conducted so close to 


the wall so the important part of the flow was fully mapped. 
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2. The Procedure 

The probes were aligned with the flow and mounted in the specially fabricated 
holder and attached to the traverse mechanism as shown in Figure 22. The traverse 
mechanism was positioned at the required points and data was recorded. The sampling 
rate was 40 KHz and 10,000 samples were taken. The data were reduced with Dantec 


Streamware software. The data taken is included in the appendices. 





Figure 22. _ Anemometer Probes 


3s Results at 2 m/s 

The entire plane was sampled at tunnel speeds of 2, 3, and 5 m/s with the 
complete results shown in Appendices D, E, and F. Figure 24 shows the turbulence 
intensity distribution at a free stream velocity of 2 m/s. Turbulence intensity ranges from 
575 % to .83 % at this speed with an average of .712 %. No particular pattern or trend is 


evident in the data. 
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Figure 23. Normalized Velocity at 2 m/s 


Another useful performance measure is the variation in velocity across the test 
section. In a perfect flow there would be no variation and no turbulence. Values of 
velocity variation across good quality wind tunnel test sections are often quoted in the .2- 
.3% range [REF 1]. The velocity variation for this facility was computed by calculating 
the percent difference between the velocity at a point and the mean velocity. The 
absolute value was then taken and these values averaged. At 2 m/s the average velocity 
variation was .251 %, well within quoted values for quality wind tunnels. Figure 23 
shows the normalized velocity distribution. Considering the scale of the axis, the data 
appear to be evenly scattered with extreme values near the top, bottom, and sides of the 


test section. 
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Figure 24. Turbulence Intensity at 2 m/s 


4. Results at 3 m/s 

At a tunnel velocity of 3 m/s, the turbulence intensity decreased to an average of 
.569% with the distribution shown in Figure 26. This downward trend in turbulence with 
an increase in velocity has been shown in past experiments. The normalized velocity 


distribution is shown in Figure 25 while the average variation at this speed was .165%. 
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Figure 25. Normalized Velocity at 3 m/s 





Turbulence Intensity at 3 m/s 


0.500 









Turbulence 
Intensity (%) 4 agg | 


0.300 — 
0.200 — 


0.100 





0.000 





w 0.800-0.900 
1 0.700-0.800 
 0.600-0.700 
@ 0.500-0.600 Y Position ( 

w 0.400-0.500 “9° = 
110.300-0.400 “8 6 

1 0.200-0.300 
 0.100-0.200 
= 0.000-0.100 





Z Position 














Figure 26. Turbulence Intensity at 3 m/s 
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5. Results at 5 m/s 


At 5 m/s the turbulence intensity crept lower to .511% and is shown in Figure 28. 


Figure 27 shows the normalized velocities with an average variation of .125%. 
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Figure 27. Normalized Velocity at 5 m/s 
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Figure 28. Turbulence Intensity at 5 m/s 


C. MODEL WAKE 


1. Flapping Wing Wake 

With a reasonable flow established in the wind tunnel, the next step was to 
analyze the MAV model wake in the facility. Initially, the model was mounted vertically 
without the main wing to see the unaltered wake of the flapping wings as shown in Figure 


29. All of the wakes shown were measured with the hotwire and are thus time averaged. 
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Figure 29. Vertical Model Without Main Wing 


The model was run at 20 Hz and 30 Hz for comparison. The hotwire was then 
traversed horizontally across the wake 10 inches on either side of the model near the 
centerline. The normalized velocity profile is shown in Figure 30 with a polynomial 
curve fit in bold. The curves are quite symmetric; especially when the polynomial fit is 
considered. This graphic clearly shows the effect of the individual wings. The velocity 
deficit between the wings at 3 m/s may indicate that the MAV must flap at a higher 
frequency to travel at 3 m/s. Alternatively, the mean separation could be lengthened to 
eliminate the deficit. Figure 31 shows the wake when flapping at 30 Hz, still indicating a 
slight velocity deficit in the center at 3 m/s. The wake appears to be somewhat narrower 
when operating at lower frequencies. Overall, the peak velocities are higher when 
flapping at 30 Hz. This illustrates the conclusion that flapping faster creates more thrust. 
In fact, thrust should increase with the square of frequency, according to the linear theory 


of waves and instabilities. 
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Figure 30. Flapping Wing Wake at 20 Hz 
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Figure 31. Flapping Wing Wake at 30 Hz 


2; Wake Behind Complete Model 
Next, the main wings were added. The model remained mounted vertically as 
shown in Figure 32. The wake was measured at a plane 19 inches from the leading edge 


of the upper flapping wing. This shifted the wake to the left side of the tunnel (positive y 
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values) in accordance with the 15° angle of attack (AOA) as shown in Figure 33. This 
AOA is typical in flight and the model was mounted in this way for a more accurate 
representation of flight. Figure 34 shows the wake shift when flapping at 30 Hz, again 
near the centerline of the model. The 3 m/s velocity deficit appears to grow when 


flapping faster. 





Figure 32. Vertical Model With Main Wings 
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Wake Behind Model at 20 Hz and 15 AOA 
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Figure 33. Wake Behind Complete Model at 20 Hz 








Wake Behind Model at 30 Hz and 15 AOA 
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Figure 34. Wake Behind Complete Model at 30 Hz 
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Figure 35. Complete Horizontal Model in the Test Section 


3. Full Wake Characterization 

Finally, the model was mounted horizontally and the complete plane of the wake 
was analyzed. The MAV was run at 30 Hz with freestream velocities of 2 m/s and 3 m/s 
and the wake was analyzed 19 inches downstream of the leading edge of the flapping 
wings. The wake was sampled in elevation from -5 to 5 inches in | inch increments. 
Each elevation was sampled in 31 points along the y axis, from -7.5 to 7.5 in .5 inch 
increments. The results from the 2 m/s run are displayed in Figure 36, with the 3 m/s 
data in Figure 37. While the flapping wings are 9.75 inches across, the wake assumes a 
different shape. At freestream velocity of 2 m/s, the bulk of the wake is concentrated in 
an area 3-5 inches wide but nearly the entire height of the test section. A similar trend is 
seen at 3 m/s with less intensity. While the flapping wings have a mean separation of 3 
inches, the wake has some level of intensity 10 inches in height. The wake appears to 
shift from short and wide coming off of the flapping wings to tall and thin a short 
distance downstream. Figure 38 illustrates the vortex structure of the wake. 
Considerable entrainment of the flow begins to occur around the leading edge of the 
flapping wings. At flight frequencies, the vortices do not change sign although they 
fluctuate in strength [REF15]. The top view of figure 38 was seen in flow visualization 
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in figure 42. The wake is drawn centerline, bringing the y axis size of the wake down. 


The aft and side views of figure 38 show the vortex structure driving the wake larger in 


the z direction as seen by the hotwire measurements. 
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Figure 37. Model Wake at 3 m/s 
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Path of vortex cores 





Side View 





Figure 38. Sketch of Vortex Structures [From REF 15] 


42 


D. FLOW VISUALIZATION 


To further investigate the flow around the model, flow visualization techniques 
were used with a smoke wire. A .075 mm tungsten wire was obtained and knots of 
copper wire were wrapped around it and soldered in place at 1 cm intervals. The ends 
were crimped in a piece of brass tubing and inserted through a special fitting. The fitting 
surrounded the wire and fed Roscoe fog fluid around the wire. The wire was strung 
vertically in the tunnel and heated with a variac. The fluid was dripped over the wire and 
pooled around the copper knots. As the wire heated, the fluid would produce smoke that 
would flow through the tunnel as streaklines. These lines lasted for several seconds. 
Experimentation revealed the voltage setting, amount of fluid on the wire, and tunnel 
speed all to be critical to successful flow visualization. The best variac setting appeared 
to be about 12 volts, and the best lines were produced after allowing the wire to cool 
briefly before starting the variac. The tunnel was operated at 2 m/s and the model was 
placed in the test section in various orientations and positions. The action was captured 


by a3 CCD Sony digital video recorder model number DCR VX-1000. 


Figure 39 shows the main wing in stall with the model stationary. The stall region 
is large and encompasses the entire upper flapping wing. Not only is this a major 
problem concerning lift, the flapping wing would not be effective in this stalled flow. 
Figure 40 shows the model operating. The flow is reattached by the entrainment that 


occurs from the flapping motion. Lift and propulsion are restored. 
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Figure 39. Main Wing Stalled When Not Flapping 





. 
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Figure 40. Flow Reattaches When Flapping 


The wingtip vortex is illustrated in Figure 41 below. Low aspect ratio wings, like 
the ones on the NPS MAV, have demonstrated an ability to resist stall. This is because 
wing tip vortices reattach the flow along the outer surfaces, restoring a significant amount 


of lift. 
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Figure 41. Wingtip Vortex 


As shown in Figures 36 and 37, the MAV wake is more intense toward the 
centerline of the craft. In the flow visualization below, the streaklines clearly move in 
this direction, creating the higher intensity wake in the vertical plane. This is also seen in 


the top view of the vortex structure of figure 38. 





Figure 42. Wake Drawn Toward Centerline 
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Vv. CONCLUSIONS 


A. TUNNEL FLOW VERIFICATION 


Ideally, a wind tunnel would produce a flow outside the boundary layer such that 
the velocity passing through each plane normal to the centerline would be uniformly 
parallel to the centerline and would have no variation with time [REF 1]. While this goal 
is not achievable in reality, real facilities need to generate flow of sufficient quality to 
fulfill their intended purpose. Values for velocity variation across the test section of high 
quality wind tunnels are often quoted in the .2-.3% range. This gives a dynamic pressure 
variation of .4 -.6% [REF 1]. The NPS MAV wind tunnel velocity variation ranged from 
125% at 5 m/s to .251% at 2 m/s. This appears to be satisfactory for MAV testing. 


B. MODEL WAKE 


The model wake shifts to a vertical orientation some distance behind it. This 


information was supported by flow visualization. 


The wingtip vortices were shown reattaching the flow to the main wings at the 
outer surfaces. This supports the MAVs ability to fly at low speeds with significant 


angles of attack. 


The action of the flapping wings acts to entrain the flow and reattach it to the 


main wings. The MAV is virtually stall proof in flight. 


C. RECOMMENDATIONS FOR FURTHER IMPROVEMENT 


Further improvements could be seen in the flow by extending the intake and 
adding more layers of screen. Some low turbulence tunnels have as many as eleven 
screen sections. In order to accomplish this addition the intake would have to be 


removed. The same mounting design could be utilized with additional hardware. The 
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frame section that supports the intake would have to be lengthened and a larger support 


fabricated to fit under the longer intake section. 


As vibration is a possible source of measurement error, the tunnel vibration could 
be reduced by remounting the motor and fan with dampers. This would limit the 
vibration transferred to the frame. Additionally, the test section exit could be replaced 


with a rubber unit that would limit the structural vibration. 
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APPENDIX A. PITOT TUBE CALIBRATION 


Voltage voltage change Tunnel speed 
(mv) (mv) (m/s) 








3.7 3.7 
8.5 48 
16 7.5 
25 9 
37 12 
47 10 
63 16 
82 19 
22 
44 
56 
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APPENDIX B. HOTWIRE CALIBRATION 


P (kpa) 


E1 corr 


U1 calc 


Fourth 
order 
polynomial 





volts 


m/s 





101.347 


1.908 


1.016 


cO 


1.977103 





101.347 


1.975 


1.32 


c1 


-8.082218 





101.344 


2.043 


1.677 


c2 


9.588953 





101.341 


2.121 


2.168 


c3 


-5.088278 





101.347 


2.208 


2.818 


c4 


1.123862 





101.352 


2.296 


3.617 





101.358 


2.393 


4.678 





101.344 


2.499 


6.087 





101.347 


2.595 


7.628 














101.347 





2.715 
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y 
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APPENDIX C. FLOW VERIFICATION AT 2 M/S 


U 
mean 
1.224 
1.511 
1.771 
2.044 
2.039 
2.047 
2.053 
2.057 
2.054 
2.044 
2.037 
2.068 
2.064 
2.062 
2.059 
2.054 
2.049 
2.051 
2.050 
2.047 
2.039 
2.051 
2.046 
2.044 
2.042 
2.051 
2.046 
2.057 
2.054 
2.044 
2.053 
2.051 


U/Uce 
0.597 
0.737 
0.864 
0.997 
0.994 
0.998 
1.001 
1.003 
1.002 
0.997 
0.993 
1.008 
1.006 
1.006 
1.004 
1.002 
0.999 
1.000 
1.000 
0.998 
0.994 
1.000 
0.998 
0.997 
0.996 
1.000 
0.998 
1.003 
1.002 
0.997 
1.001 


T% 
1.316 
1.125 
0.960 
0.864 
0.854 
0.825 
0.781 
0.765 
0.779 
0.811 
0.757 
0.778 
0.738 
0.730 
0.741 
0.761 
0.784 
0.746 
0.789 
0.753 
0.741 
0.755 
0.793 
0.785 
0.742 
0.763 
0.802 
0.769 
0.784 
0.811 
0.774 
AVG 


U var 
% 
40.313 
26.318 
13.639 
0.327 
0.570 
0.180 
0.112 
0.307 
0.161 
0.327 
0.668 
0.844 
0.649 
0.551 
0.405 
0.161 
0.083 
0.015 
0.034 
0.180 
0.570 
0.015 
0.229 
0.327 
0.424 
0.015 
0.229 
0.307 
0.161 
0.327 
0.112 
0.301 
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U/Us 
0.567 
0.758 
0.855 
0.947 
1.007 
1.006 
1.010 
1.010 
1.007 
1.006 
1.006 
1.007 
1.006 
1.003 
1.003 
1.003 
1.003 
1.001 
1.000 
1.001 
0.999 
0.998 
0.997 
0.998 
0.995 
0.995 
0.990 
0.990 
0.995 
0.997 
0.999 


T% 
1.411 
1.161 
1.027 
0.865 
0.819 
0.835 
0.866 
0.770 
0.682 
0.741 
0.779 
0.722 
0.763 
0.761 
0.672 
0.731 
0.601 
0.660 
0.752 
0.656 
0.721 
0.615 
0.661 
0.646 
0.641 
0.648 
0.725 
0.736 
0.620 
0.691 
0.765 
AVG 


U var 
% 
43.298 
24.235 
14.484 
5.318 
0.678 
0.630 
1.020 
0.971 
0.727 
0.581 
0.630 
0.678 
0.581 
0.288 
0.337 
0.337 
0.337 
0.093 
0.004 
0.093 
0.102 
0.199 
0.297 
0.199 
0.492 
0.492 
0.979 
0.979 
0.541 
0.297 
0.102 
0.393 


y 
11.875 


11.813 
11.750 
11.688 
11.625 
11.500 
11.375 
11.250 
11.125 
11.000 
10.000 
9.000 
8.000 
7.000 
6.000 
5.000 
4.000 
3.000 
2.000 
1.000 
0.000 
-1.000 
-2.000 
-3.000 
-4.000 
-5.000 
-6.000 
-7.000 
-8.000 
-9.000 
-10.000 
Uo= 


U 
mean 
1.202 
1.405 
1.703 
1.807 
2.032 
2.069 
2.071 
2.067 
2.069 
2.068 
2.066 
2.064 
2.064 
2.060 
2.059 
2.057 
2.055 
2.058 
2.055 
2.052 
2.055 
2.054 
2.053 
2.054 
2.058 
2.050 
2.050 
2.052 
2.045 
2.043 
2.041 
2.056 
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U/Uce 
0.579 
0.721 
0.835 
0.940 
0.976 
0.997 
1.004 
1.003 
1.004 
1.003 
1.004 
1.002 
1.000 
1.000 
1.001 
0.999 
0.998 
0.996 
0.999 
0.997 
0.997 
0.998 
0.999 
0.996 
0.999 
0.999 
1.000 
1.002 
1.002 
1.003 
1.002 


T% 
1.111 
1.013 
0.816 
0.760 
0.828 
0.674 
0.729 
0.743 
0.746 
0.748 
0.603 
0.756 
0.756 
0.661 
0.612 
0.641 
0.645 
0.662 
0.613 
0.661 
0.652 
0.711 
0.621 
0.752 
0.655 
0.732 
0.635 
0.622 
0.611 
0.603 
0.610 
AVG 


U var 
% 
42.051 
27.879 
16.522 
5.991 
2.351 
0.312 
0.416 
0.319 
0.367 
0.319 
0.367 
0.173 
0.021 
0.021 
0.076 
0.118 
0.215 
0.361 
0.070 
0.264 
0.264 
0.215 
0.118 
0.361 
0.118 
0.070 
0.027 
0.173 
0.173 
0.319 
0.222 
0.185 


Z 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 


y 
11.875 


11.813 
11.750 
11.688 
11.625 
11.500 
11.375 
11.250 
11.125 
11.000 
10.000 
9.000 
8.000 
7.000 
6.000 
5.000 
4.000 
3.000 
2.000 
1.000 
0.000 
-1.000 
-2.000 
-3.000 
-4.000 
-5.000 
-6.000 
-7.000 
-8.000 
-9.000 
-10.000 
Uo= 


U/Use 
0.539 
0.710 
0.846 
0.925 
0.976 
1.006 
1.007 
1.005 
1.008 
1.006 
1.005 
1.004 
1.002 
1.002 
1.000 
0.998 
0.998 
0.997 
0.997 
0.998 
0.997 
0.998 
0.998 
0.997 
0.999 
0.999 
1.001 
0.998 
0.998 
0.997 
0.998 


T% 
1.019 
0.926 
0.839 
0.761 
0.706 
0.680 
0.741 
0.704 
0.710 
0.855 
0.740 
0.700 
0.762 
0.695 
0.848 
0.800 
0.734 
0.815 
0.668 
0.726 
0.653 
0.713 
0.637 
0.701 
0.752 
0.786 
0.633 
0.651 
0.701 
0.734 
0.746 
AVG 


U var 
% 
46.052 
28.980 
15.371 
7.469 
2.396 
0.579 
0.677 
0.530 
0.774 
0.579 
0.530 
0.433 
0.238 
0.189 
0.042 
0.153 
0.201 
0.250 
0.250 
0.201 
0.250 
0.201 
0.153 
0.299 
0.104 
0.055 
0.091 
0.153 
0.201 
0.250 
0.153 
0.226 


55 


Z 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 


y 
11.875 


11.813 
11.750 
11.688 
11.625 
11.500 
11.375 
11.250 
11.125 
11.000 
10.000 
9.000 
8.000 
7.000 
6.000 
5.000 
4.000 
3.000 
2.000 
1.000 
0.000 
-1.000 
-2.000 
-3.000 
-4.000 
-5.000 
-6.000 
-7.000 
-8.000 
-9.000 
-10.000 
Uo= 


U/Use 
0.543 
0.694 
0.851 
0.922 
0.982 
1.002 
0.994 
0.997 
0.995 
1.005 
1.003 
1.002 
1.000 
1.001 
1.002 
0.999 
0.998 
1.001 
1.000 
0.999 
0.999 
1.000 
0.999 
0.997 
1.000 
1.001 
0.999 
1.001 
0.998 
0.999 
0.998 


T% 
1.019 
0.926 
0.839 
0.740 
0.695 
0.682 
0.736 
0.733 
0.735 
0.873 
0.729 
0.682 
0.780 
0.682 
0.827 
0.781 
0.733 
0.828 
0.683 
0.732 
0.634 
0.731 
0.635 
0.685 
0.780 
0.779 
0.634 
0.649 
0.701 
0.735 
0.746 
AVG 


U var 
% 
45.667 
30.596 
14.891 
7.770 
1.819 
0.180 
0.551 
0.259 
0.454 
0.522 
0.327 
0.180 
0.034 
0.131 
0.229 
0.064 
0.161 
0.083 
0.034 
0.064 
0.064 
0.034 
0.112 
0.307 
0.015 
0.131 
0.064 
0.083 
0.161 
0.112 
0.210 
0.142 


z 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 


y 
11.875 


11.813 
11.750 
11.688 
11.625 
11.500 
11.375 
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11.125 
11.000 
10.000 
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8.000 
7.000 
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4.000 
3.000 
2.000 
1.000 
0.000 
-1.000 
-2.000 
-3.000 
-4.000 
-5.000 
-6.000 
-7.000 
-8.000 
-9.000 
-10.000 
Uo= 


U mean 
1.122 
1.512 
1.743 
1.929 
2.010 
2.052 
2.061 
2.042 
2.056 
2.058 
2.047 
2.049 
2.055 
2.059 
2.044 
2.043 
2.051 
2.059 
2.052 
2.049 
2.056 
2.053 
2.062 
2.041 
2.045 
2.040 
2.047 
2.054 
2.051 
2.045 
2.047 
2.051 


Overall velocity variation at 2 m/s 
(%) 


Average turbulence intensity at 2 
m/s (%) 
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U/Use 
0.547 
0.737 
0.850 
0.941 
0.980 
1.001 
1.005 
0.996 
1.003 
1.004 
0.998 
0.999 
1.002 
1.004 
0.997 
0.996 
1.000 
1.004 
1.001 
0.999 
1.003 
1.001 
1.006 
0.995 
0.997 
0.995 
0.998 
1.002 
1.000 
0.997 
0.998 


0.980 
0.860 
0.861 
0.726 
0.697 
0.634 
0.728 
0.735 
0.681 
0.715 
0.733 
0.683 
0.779 
0.680 
0.732 
0.783 
0.731 
0.777 
0.682 
0.732 
0.632 
0.731 
0.630 
0.735 
0.782 
0.784 
0.635 
0.779 
0.731 
0.782 
0.830 
AVG 


APPENDIX D. FLOW VERIFICATION AT 3 M/S 


U/Use 
0.833 
0.900 
0.946 
0.968 
0.964 
0.994 
0.997 
1.002 
1.002 
1.003 
1.001 
1.005 
1.001 
1.000 
1.001 
1.000 
1.000 
1.000 
0.999 
1.001 
0.999 
0.999 
0.999 
0.998 
0.998 
1.000 
0.999 
0.999 
0.999 
1.000 
0.999 


T% 
11.035 
8.528 
5.938 
4.424 
5.411 
2.358 
1.932 
1.416 
1.194 
1.082 
0.634 
0.600 
0.571 
0.607 
0.563 
0.496 
0.564 
0.562 
0.556 
0.551 
0.581 
0.587 
0.506 
0.632 
0.536 
0.549 
0.523 
0.513 
0.566 
0.564 
0.573 
AVG 


U var 
% 
16.742 
9.963 
5.411 
3.248 
3.603 
0.633 
0.278 
0.239 
0.207 
0.336 
0.077 
0.465 
0.142 
0.013 
0.077 
0.019 
0.019 
0.013 
0.084 
0.077 
0.084 
0.116 
0.149 
0.181 
0.213 
0.045 
0.149 
0.149 
0.084 
0.019 
0.149 
0.121 
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U/Us 
0.878 
0.952 
1.002 
0.992 
1.004 
1.005 
1.004 
1.005 
1.006 
1.004 
1.002 
1.003 
1.002 
1.000 
0.999 
0.997 
0.998 
0.998 
0.997 
0.998 
0.999 
0.997 
0.999 
0.999 
0.999 
0.998 
0.998 
0.999 
0.999 
0.999 
0.999 


T% 
10.026 
5.544 
1.297 
3.761 
0.907 
0.782 
0.903 
0.701 
0.658 
0.560 
0.552 
0.589 
0.532 
0.579 
0.543 
0.495 
0.530 
0.528 
0.562 
0.528 
0.518 
0.511 
0.600 
0.505 
0.615 
0.512 
0.535 
0.573 
0.547 
0.568 
0.556 
AVG 


U var 
% 
12.164 
4.834 
0.172 
0.829 
0.398 
0.462 
0.398 
0.462 
0.559 
0.430 
0.236 
0.333 
0.204 
0.043 
0.119 
0.313 
0.183 
0.183 
0.280 
0.183 
0.119 
0.280 
0.119 
0.119 
0.054 
0.183 
0.216 
0.119 
0.087 
0.054 
0.054 
0.178 


U/Use 
0.848 
0.896 
0.960 
0.998 
0.999 
0.991 
1.004 
1.004 
1.003 
1.002 
1.003 
1.001 
1.002 
0.999 
0.999 
1.000 
1.000 
0.999 
0.998 
0.998 
0.998 
0.999 
0.998 
0.999 
1.000 
0.999 
0.999 
0.999 
0.999 
0.999 
0.999 
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U/Us 
0.797 
0.834 
0.907 
0.898 
0.944 
0.959 
0.995 
1.000 
1.008 
1.007 
1.006 
1.005 
1.005 
1.002 
1.002 
1.002 
1.001 
1.000 
0.999 
0.998 
0.998 
0.998 
0.998 
0.998 
0.997 
0.998 
0.997 
0.998 
0.997 
0.994 
0.994 


T% 
13.109 
9.389 
6.041 
6.862 
6.001 
5.507 
2.117 
2.152 
0.800 
0.642 
0.614 
0.585 
0.599 
0.529 
0.570 
0.513 
0.542 
0.611 
0.581 
0.641 
0.557 
0.579 
0.585 
0.597 
0.531 
0.539 
0.527 
0.553 
0.475 
0.599 
0.549 
AVG 


U var 
% 
20.307 
16.614 
9.292 
10.199 
5.632 
4.077 
0.546 
0.038 
0.750 
0.653 
0.621 
0.524 
0.491 
0.232 
0.200 
0.200 
0.102 
0.038 
0.060 
0.157 
0.157 
0.157 
0.189 
0.189 
0.286 
0.189 
0.254 
0.189 
0.286 
0.578 
0.610 
0.289 


Z 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 


y 
11.875 


11.813 
11.750 
11.688 
11.625 
11.500 
11.375 
11.250 
11.125 
11.000 
10.000 
9.000 
8.000 
7.000 
6.000 
5.000 
4.000 
3.000 
2.000 
1.000 
0.000 
-1.000 
-2.000 
-3.000 
-4.000 
-5.000 
-6.000 
-7.000 
-8.000 
-9.000 
-10.000 
Uo= 


U/Use 
0.753 
0.839 
0.861 
0.900 
0.900 
0.957 
0.987 
1.001 
1.004 
1.003 
1.002 
1.002 
1.000 
1.001 
1.000 
1.001 
1.002 
1.001 
1.003 
0.998 
1.000 
1.001 
1.000 
1.001 
0.998 
1.000 
1.001 
0.999 
1.000 
0.999 
0.999 


T% 
12.631 
7.396 
8.497 
6.360 
5.505 
4.361 
2.302 
1.096 
0.848 
0.791 
0.566 
0.580 
0.576 
0.553 
0.544 
0.519 
0.548 
0.556 
0.653 
0.622 
0.533 
0.540 
0.566 
0.515 
0.569 
0.609 
0.523 
0.629 
0.557 
0.600 
0.627 
AVG 


U var 
% 
24.652 
16.123 
13.934 
10.007 
9.975 
4.343 
1.317 
0.099 
0.357 
0.260 
0.164 
0.228 
0.035 
0.131 
0.030 
0.099 
0.164 
0.067 
0.260 
0.158 
0.030 
0.067 
0.030 
0.067 
0.158 
0.035 
0.067 
0.126 
0.030 
0.094 
0.126 
0.110 
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Z 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 


y 
11.875 


11.813 
11.750 
11.688 
11.625 
11.500 
11.375 
11.250 
11.125 
11.000 
10.000 
9.000 
8.000 
7.000 
6.000 
5.000 
4.000 
3.000 
2.000 
1.000 
0.000 
-1.000 
-2.000 
-3.000 
-4.000 
-5.000 
-6.000 
-7.000 
-8.000 
-9.000 
-10.000 
Uo= 


U/Us 
0.855 
0.850 
0.892 
0.954 
0.973 
0.999 
1.002 
1.007 
1.004 
1.004 
1.002 
1.002 
1.001 
1.002 
1.001 
1.001 
1.001 
1.001 
1.001 
1.001 
1.000 
1.000 
1.001 
0.997 
0.997 
0.997 
0.997 
0.998 
0.997 
0.994 
0.992 


T% 
9.322 
7.983 
6.502 
4.936 
3.703 
1.831 
1.396 
0.823 
0.563 
0.601 
0.572 
0.526 
0.561 
0.544 
0.536 
0.532 
0.582 
0.548 
0.508 
0.550 
0.533 
0.604 
0.556 
0.579 
0.521 
0.600 
0.589 
0.532 
0.558 
0.537 
0.554 
AVG 


U var 
% 
14.474 
15.021 
10.806 
4.564 
2.698 
0.091 
0.230 
0.713 
0.423 
0.391 
0.230 
0.166 
0.070 
0.198 
0.102 
0.134 
0.102 
0.134 
0.070 
0.070 
0.005 
0.005 
0.070 
0.349 
0.284 
0.317 
0.284 
0.188 
0.252 
0.638 
0.799 
0.221 


z 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 


y 
11.875 


11.813 
11.750 
11.688 
11.625 
11.500 
11.375 
11.250 
11.125 
11.000 
10.000 
9.000 
8.000 
7.000 
6.000 
5.000 
4.000 
3.000 
2.000 
1.000 
0.000 
-1.000 
-2.000 
-3.000 
-4.000 
-5.000 
-6.000 
-7.000 
-8.000 
-9.000 
-10.000 
Uo= 


U mean 
2.284 
2.317 
2.413 
2.463 
2.618 
2.704 
2.819 
2.893 
2.852 
2.969 
3.064 
3.069 
3.065 
3.063 
3.068 
3.062 
3.063 
3.068 
3.064 
3.059 
3.061 
3.056 
3.059 
3.064 
3.065 
3.063 
3.069 
3.056 
3.054 
3.058 
3.060 
3.062 


Overall velocity variation at 3 m/s 
(%) 


Average turbulence intensity at 3 
m/s (%) 
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U/Use 
0.746 
0.757 
0.788 
0.804 
0.855 
0.883 
0.921 
0.945 
0.931 
0.970 
1.001 
1.002 
1.001 
1.000 
1.002 
1.000 
1.000 
1.002 
1.001 
0.999 
1.000 
0.998 
0.999 
1.001 
1.001 
1.000 
1.002 
0.998 
0.997 
0.999 
0.999 


T% 
9.466 
10.734 
9.987 
8.254 
7.596 
6.997 
5.770 
5.859 
6.649 
4.363 
0.875 
0.618 
0.559 
0.532 
0.565 
0.572 
0.622 
0.546 
0.596 
0.631 
0.599 
0.634 
0.608 
0.570 
0.581 
0.659 
0.580 
0.648 
0.614 
0.579 
0.574 
AVG 


U var 
% 
25.418 
24.340 
21.205 
19.572 
14.511 
11.703 
7.947 
5.531 
6.870 
3.049 
0.053 
0.216 
0.086 
0.020 
0.183 
0.012 
0.020 
0.183 
0.053 
0.110 
0.045 
0.208 
0.110 
0.053 
0.086 
0.020 
0.216 
0.208 
0.274 
0.143 
0.078 
0.113 


0.165 


0.569 


APPENDIX E. FLOW VERIFICATION AT 5 M/S 


U/Use 
0.857 
0.888 
0.931 
0.941 
0.968 
0.991 
0.998 
0.997 
0.997 
1.000 
0.996 
0.998 
0.998 
0.996 
0.997 
1.001 
1.001 
1.002 
1.001 
1.001 
1.002 
1.002 
1.001 
1.002 
1.001 
1.000 
1.000 
1.001 
1.000 
1.001 
1.002 
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U/Use 
0.892 
0.909 
0.942 
0.966 
0.982 
1.003 
1.003 
1.002 
1.004 
1.005 
1.004 
1.003 
1.003 
1.001 
1.000 
1.000 
0.998 
1.000 
0.999 
0.999 
0.998 
0.998 
0.998 
0.999 
0.997 
0.997 
1.000 
0.999 
0.998 
0.999 
0.998 


T% 
6.457 
5.801 
4.757 
4.187 
3.422 
1.007 
0.880 
0.714 
0.603 
0.561 
0.502 
0.564 
0.466 
0.515 
0.551 
0.510 
0.518 
0.514 
0.563 
0.554 
0.578 
0.496 
0.540 
0.529 
0.516 
0.515 
0.536 
0.475 
0.511 
0.508 
0.587 
AVG 


U var 
% 
10.792 
9.130 
5.844 
3.400 
1.757 
0.315 
0.276 
0.218 
0.394 
0.491 
0.413 
0.276 
0.257 
0.100 
0.017 
0.002 
0.174 
0.037 
0.095 
0.115 
0.174 
0.174 
0.174 
0.134 
0.271 
0.252 
0.037 
0.076 
0.174 
0.056 
0.193 
0.168 


y 
11.875 


11.813 
11.750 
11.688 
11.625 
11.500 
11.375 
11.250 
11.125 
11.000 
10.000 
9.000 
8.000 
7.000 
6.000 
5.000 
4.000 
3.000 
2.000 
1.000 
0.000 
-1.000 
-2.000 
-3.000 
-4.000 
-5.000 
-6.000 
-7.000 
-8.000 
-9.000 
-10.000 
Uo= 


U 
mean 
4.340 
4.533 
4.719 
4.903 
4.992 
5.111 
5.146 
5.145 
5.146 
5.148 
5.146 
5.144 
5.141 
5.133 
5.133 
5.131 
5.136 
5.130 
5.122 
5.129 
5.121 
5.126 
5.119 
5.120 
5.120 
5.120 
5.125 
5.128 
5.116 
5.124 
5.128 
5.130 


U/Use 
0.846 
0.884 
0.920 
0.956 
0.973 
0.996 
1.003 
1.003 
1.003 
1.003 
1.003 
1.003 
1.002 
1.000 
1.000 
1.000 
1.001 
1.000 
0.998 
1.000 
0.998 
0.999 
0.998 
0.998 
0.998 
0.998 
0.999 
1.000 
0.997 
0.999 
1.000 
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U/Use 
0.832 
0.884 
0.890 
0.937 
0.940 
0.981 
0.997 
0.999 
1.001 
1.002 
0.999 
1.001 
1.002 
1.000 
1.000 
1.000 
1.000 
0.999 
0.999 
1.000 
1.000 
1.000 
1.000 
1.001 
1.000 
0.999 
0.999 
1.000 
1.000 
1.000 
1.000 


T% 
8.802 
7.019 
6.536 
4.842 
5.334 
2.656 
1.059 
0.672 
0.675 
0.557 
0.471 
0.542 
0.525 
0.539 
0.476 
0.500 
0.529 
0.491 
0.473 
0.528 
0.475 
0.533 
0.562 
0.615 
0.507 
0.483 
0.438 
0.461 
0.463 
0.543 
0.483 
AVG 


U var 
% 
16.848 
11.600 
11.039 
6.256 
5.984 
1.918 
0.291 
0.098 
0.096 
0.154 
0.059 
0.057 
0.154 
0.040 
0.040 
0.038 
0.019 
0.078 
0.098 
0.038 
0.040 
0.038 
0.019 
0.077 
0.001 
0.078 
0.117 
0.020 
0.020 
0.040 
0.038 
0.057 


Z 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 


y 
11.875 


11.813 
11.750 
11.688 
11.625 
11.500 
11.375 
11.250 
11.125 
11.000 
10.000 
9.000 
8.000 
7.000 
6.000 
5.000 
4.000 
3.000 
2.000 
1.000 
0.000 
-1.000 
-2.000 
-3.000 
-4.000 
-5.000 
-6.000 
-7.000 
-8.000 
-9.000 
-10.000 
Uo= 


U 
mean 
4.182 
4.284 
4.684 
4.830 
4.876 
5.121 
5.153 
5.167 
5.158 
5.162 
5.163 
5.157 
5.149 
5.149 
5.148 
5.149 
5.142 
5.146 
5.149 
5.152 
5.149 
5.156 
5.154 
5.154 
5.150 
5.153 
5.153 
5.156 
5.159 
5.154 
5.158 
5.154 


U/Use 
0.811 
0.831 
0.909 
0.937 
0.946 
0.994 
1.000 
1.003 
1.001 
1.002 
1.002 
1.001 
0.999 
0.999 
0.999 
0.999 
0.998 
0.999 
0.999 
1.000 
0.999 
1.000 
1.000 
1.000 
0.999 
1.000 
1.000 
1.000 
1.001 
1.000 
1.001 


T% 
8.625 
7.406 
6.042 
5.217 
4.511 
1.907 
1.066 
0.634 
0.626 
0.508 
0.490 
0.503 
0.510 
0.490 
0.523 
0.462 
0.453 
0.474 
0.535 
0.456 
0.518 
0.502 
0.462 
0.461 
0.512 
0.477 
0.565 
0.568 
0.470 
0.509 
0.461 
AVG 


U var 
% 
18.853 
16.874 
9.113 
6.280 
5.387 
0.633 
0.012 
0.260 
0.085 
0.163 
0.182 
0.065 
0.090 
0.090 
0.109 
0.090 
0.226 
0.148 
0.090 
0.032 
0.090 
0.046 
0.007 
0.007 
0.070 
0.012 
0.012 
0.046 
0.104 
0.007 
0.085 
0.080 
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Z 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 


y 
11.875 


11.813 
11.750 
11.688 
11.625 
11.500 
11.375 
11.250 
11.125 
11.000 
10.000 
9.000 
8.000 
7.000 
6.000 
5.000 
4.000 
3.000 
2.000 
1.000 
0.000 
-1.000 
-2.000 
-3.000 
-4.000 
-5.000 
-6.000 
-7.000 
-8.000 
-9.000 
-10.000 
Uo= 


U/Uce 
0.792 
0.844 
0.879 
0.905 
0.930 
0.971 
0.998 
1.001 
1.001 
1.003 
1.001 
1.000 
1.001 
0.999 
1.000 
0.999 
1.000 
0.999 
0.999 
1.000 
1.000 
0.999 
1.000 
0.999 
1.000 
1.000 
1.000 
1.000 
0.998 
1.001 
1.000 


T% 
8.130 
8.658 
6.125 
6.350 
5.738 
3.575 
1.557 
0.878 
0.630 
0.555 
0.536 
0.465 
0.515 
0.480 
0.480 
0.476 
0.480 
0.522 
0.430 
0.473 
0.438 
0.480 
0.581 
0.513 
0.521 
0.472 
0.500 
0.422 
0.471 
0.498 
0.542 
AVG 


U var 
% 
20.775 
15.626 
12.103 
9.529 
6.955 
2.948 
0.160 
0.130 
0.091 
0.285 
0.149 
0.014 
0.072 
0.064 
0.014 
0.122 
0.025 
0.083 
0.064 
0.025 
0.033 
0.064 
0.014 
0.083 
0.044 
0.025 
0.044 
0.006 
0.238 
0.091 
0.006 
0.071 


z 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 
-6.0 


y 
11.875 


11.813 
11.750 
11.688 
11.625 
11.500 
11.375 
11.250 
11.125 
11.000 
10.000 
9.000 
8.000 
7.000 
6.000 
5.000 
4.000 
3.000 
2.000 
1.000 
0.000 
-1.000 
-2.000 
-3.000 
-4.000 
-5.000 
-6.000 
-7.000 
-8.000 
-9.000 
-10.000 
Uo= 


U mean 
4.087 
4.237 
4.430 
4.593 
4.675 
4.785 
5.007 
5.159 
5.151 
5.184 
5.177 
5.179 
5.176 
5.169 
5.156 
5.163 
5.160 
5.149 
5.146 
5.152 
5.151 
5.157 
5.151 
5.153 
5.155 
5.147 
5.155 
5.145 
5.160 
5.150 
5.148 
5.158 


Overall velocity variation at 5 m/s 
(%) 


Average turbulence intensity at 5 
m/s (%) 
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U/Uce 
0.792 
0.821 
0.859 
0.890 
0.906 
0.928 
0.971 
1.000 
0.999 
1.005 
1.004 
1.004 
1.003 
1.002 
1.000 
1.001 
1.000 
0.998 
0.998 
0.999 
0.999 
1.000 
0.999 
0.999 
0.999 
0.998 
0.999 
0.997 
1.000 
0.998 
0.998 


T% 
9.185 
9.664 
7.971 
7.468 
6.892 
6.167 
4.415 
1.882 
1.911 
1.200 
0.550 
0.539 
0.485 
0.555 
0.455 
0.510 
0.548 
0.489 
0.499 
0.599 
0.526 
0.543 
0.454 
0.561 
0.589 
0.544 
0.458 
0.566 
0.462 
0.467 
0.550 
AVG 


U var 
% 
20.765 
17.856 
14.115 
10.955 
9.365 
7.232 
2.928 
0.019 
0.137 
0.503 
0.368 
0.406 
0.348 
0.212 
0.040 
0.096 
0.038 
0.175 
0.233 
0.117 
0.137 
0.020 
0.137 
0.098 
0.059 
0.214 
0.059 
0.253 
0.038 
0.156 
0.195 
0.177 


0.125 


0.511 


APPENDIX F. FLAPPING WING WAKE AT CONSTANT Z 


U/Us 
1.017 
1.018 
1.018 
1.097 
1.235 
1.276 
1.154 
1.164 
1.124 
1.185 
1.068 
1.098 
1.138 
1.263 
1.270 
1.148 
0.990 
0.984 
0.959 
1.014 


U RMS 
0.016 
0.026 
0.087 
0.210 
0.234 
0.234 
0.300 
0.269 
0.260 
0.236 
0.260 
0.324 
0.359 
0.389 
0.354 
0.402 
0.216 
0.141 
0.075 
0.047 


POSITION 
T% Zz 
0.729 -1.0 
1.186 -1.0 
4.044 -1.0 
9.051 -1.0 
8.923 -1.0 
8.651 -1.0 
12.283 -1.0 
10.919 -1.0 
10.929 -1.0 
9.383 -1.0 
11.488 -1.0 
13.948 -1.0 
14.869 -1.0 
14.523 -1.0 
13.176 -1.0 
12.520 -1.0 
10.319 -1.0 
6.769 -1.0 
3.666 -1.0 
2.195 -1.0 


65 


U/Us 
1.002 
1.001 
1.001 
1.001 
1.138 
1.045 
1.055 
1.089 
0.966 
0.905 
0.970 
1.141 
1.207 
1.108 
1.033 
1.010 
0.974 
0.902 
0.973 
0.997 


U RMS 
0.018 
0.019 
0.025 
0.050 
0.143 
0.435 
0.357 
0.340 
0.322 
0.319 
0.342 
0.259 
0.363 
0.419 
0.376 
0.256 
0.106 
0.096 
0.097 
0.051 


T% 
0.582 
0.607 
0.823 
1.615 
4.613 

14.112 
11.012 
10.180 
10.859 
11.456 
12.318 
10.385 
11.799 
12.308 
11.836 
8.267 
3.555 
3.457 
3.232 
1.658 
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APPENDIX G. FULL MODEL WAKE AT 30 HZ AND 2 M/S 


U/Us 
1.009 
1.007 
1.007 
1.007 
1.006 
1.005 
1.004 
1.003 
1.004 
1.001 
1.001 
1.011 
1.094 
1.075 
1.172 
1.224 
1.179 
1.102 
1.027 
1.012 
1.000 
0.999 
1.001 
0.999 
0.999 
0.998 
0.996 
0.995 
0.996 
0.994 
0.985 
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U/U 
1.005 
1.005 
1.003 
1.004 
1.008 
1.006 
1.003 
1.001 
1.000 
1.002 
1.004 
1.033 
1.165 
1.234 
1.277 
1.324 
1.224 
1.238 
1.072 
1.011 
1.003 
1.002 
1.001 
0.998 
1.001 
1.001 
0.999 
0.995 
0.995 
0.991 
0.986 


T% 
0.894 
0.973 
0.958 
0.812 
0.861 
0.940 
0.832 
0.910 
1.121 
1.807 
2.664 
8.218 
12.090 
12.432 
9.223 
10.678 
11.079 
9.569 
9.551 
4.463 
1.491 
1.349 
0.955 
0.695 
0.746 
0.743 
0.763 
0.772 
0.784 
0.866 
0.702 


U/U% 
1.005 
1.007 
1.008 
1.008 
1.009 
1.010 
1.009 
1.009 
1.000 
1.003 
1.003 
1.037 
1.072 
1.170 
1.374 
1.270 
1.278 
1.182 
1.032 
1.002 
0.999 
1.000 
1.000 
0.999 
1.000 
0.998 
0.999 
0.996 
0.993 
0.991 
0.987 
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U/U% 
1.006 
1.010 
1.004 
1.003 
1.010 
1.017 
1.014 
1.006 
1.000 
0.998 
0.998 
0.998 
1.006 
1.160 
1.181 
1.347 
1.296 
1.072 
1.008 
1.010 
0.997 
1.000 
1.002 
1.004 
1.000 
1.003 
0.999 
0.997 
0.994 
0.989 
0.985 


T% 
1.276 
1.016 
0.911 
1.755 
2.183 
3.316 
2.535 
1.488 
2.353 
1.997 
1.748 
2.422 
5.741 
16.450 
14.004 
12.645 
13.003 
10.830 
6.082 
6.605 
1.522 
1.272 
1.704 
1.104 
0.916 
1.554 
1.321 
1.004 
0.785 
0.690 
0.701 


U/U 
1.042 
1.058 
1.018 
1.019 
1.011 
1.010 
1.018 
1.051 
1.060 
1.081 
1.042 
0.993 
1.006 
1.000 
1.190 
1.385 
1.274 
1.067 
0.999 
1.002 
1.009 
1.020 
1.014 
1.032 
1.022 
1.007 
1.001 
1.028 
1.002 
0.997 
0.979 
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U/U% 
1.018 
1.029 
1.004 
0.975 
0.994 
0.976 
1.007 
1.051 
1.084 
1.080 
1.073 
1.085 
1.032 
1.009 
1.164 
1.316 
1.180 
1.033 
0.973 
1.015 
1.032 
1.093 
1.095 
1.043 
1.054 
1.054 
0.998 
1.008 
1.022 
0.986 
0.983 


T% 
4.540 
6.992 
6.199 
6.703 
3.720 
5.669 
4.132 
5.564 
4.559 
5.371 
7.515 
5.335 
9.574 
12.536 
13.186 
11.574 
9.801 
11.536 
7.350 
4.398 
8.192 
7427 
7.217 
8.131 
8.485 
7.978 
5.160 
4.490 
6.250 
1.206 
1.061 


U/U% 
1.010 
1.018 
1.010 
1.000 
0.994 
1.008 
1.012 
1.023 
1.036 
1.113 
1.088 
1.059 
0.996 
1.066 
1.146 
1.302 
1.237 
1.013 
1.014 
1.079 
1.077 
1.101 
1.078 
1.034 
1.057 
1.007 
1.011 
1.009 
1.000 
0.986 
0.983 


T% 
2.436 
3.418 
3.128 
3.632 
3.693 
3.043 
4.340 
3.204 
7.038 
10.507 
10.509 
8.794 
11.085 
13.873 
14.149 
13.664 
12.194 
10.530 
10.043 
10.024 
8.985 
7.945 
8.191 
5.298 
8.011 
8.612 
7.016 
4.394 
6.801 
2.126 
1.074 
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-2.0 


U/U 
1.002 
1.005 
1.005 
1.007 
1.002 
1.008 
1.002 
1.007 
1.046 
1.132 
1.226 
1.159 
1.080 
1.148 
1.279 
1.261 
1.098 
1.093 
1.060 
1.141 
1.102 
1.031 
1.044 
1.031 
1.020 
0.998 
1.021 
1.008 
0.999 
0.988 
0.984 


T% 
0.837 
0.919 
1.296 
0.916 
1.968 
2.089 
3.129 
6.143 
8.464 
9.448 
9.815 

10.561 
13.458 
11.367 
10.999 
11.560 
12.002 
9.629 
12.151 
11.267 
10.769 
8.024 
6.946 
6.287 
6.362 
6.474 
6.652 
3.669 
3.159 
1.494 
0.934 


U/Us 
1.005 
1.009 
1.011 
1.011 
1.017 
1.020 
1.013 
1.007 
1.010 
1.020 
1.159 
1.226 
1.151 
1.163 
1.187 
1.229 
1.247 
1.182 
1.170 
1.141 
1.049 
1.064 
1.023 
1.049 
1.025 
1.020 
1.011 
1.008 
0.998 
0.993 
0.987 


T% 
0.715 
0.902 
1.086 
1.119 
1.344 
1.383 
3.384 
6.436 
6.073 
7.953 
10.293 
9.178 
9.243 
10.629 
9.613 
11.628 
12.095 
9.694 
8.800 
12.434 
9.832 
8.222 
7.236 
7.976 
4.477 
4.096 
3.523 
1.934 
1.281 
0.970 
0.867 
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-4.0 


U/Us 
1.004 
1.020 
1.023 
1.027 
1.027 
1.030 
1.024 
1.013 
0.984 
1.008 
1.032 
1.093 
1.080 
1.073 
1.079 
1.090 
1.039 
1.025 
1.012 
1.025 
1.024 
1.003 
1.009 
1.045 
1.028 
1.023 
1.005 
1.007 
0.991 
0.991 
0.983 


T% 
1.226 
1.384 
1.967 
2.238 
1.035 
1.445 
2.860 
4.365 
9.445 
7.696 
6.546 
7.890 
11.087 
9.076 
11.032 
13.037 
10.511 
12.061 
9.184 
10.401 
7.059 
8.828 
7.553 
3.749 
1.940 
3.642 
1.969 
1.418 
2.473 
0.863 
1.035 


U/Us 
0.984 
0.989 
1.029 
1.034 
1.054 
1.060 
1.057 
1.050 
1.066 
1.012 
1.028 
1.049 
1.224 
1.181 
1.295 
1.305 
1.187 
1.035 
1.059 
1.040 
1.043 
0.996 
1.028 
1.049 
1.064 
1.048 
1.007 
0.993 
1.001 
0.966 
0.996 


T% 
5.122 
5.116 
5.786 
6.763 
5.739 
5.923 
5.180 
6.443 
9.042 
10.503 
10.102 
12.905 
13.142 
15.975 
11.470 
12.110 
13.422 
12.268 
12.353 
10.802 
9.264 
8.662 
5.796 
7.021 
6.412 
5.639 
6.293 
5.510 
3.583 
6.158 
2.555 


a2 


APPENDIX H. FULL MODEL WAKE AT 30 HZ AND 3 M/S 


U/Use 
1.007 
1.005 
1.012 
1.011 
1.012 
1.013 
1.012 
1.014 
1.005 
1.003 
1.003 
1.003 
1.004 
1.003 
1.007 
1.000 
1.001 
1.005 
1.003 
1.003 
1.001 
1.003 
1.002 
1.003 
1.002 
1.000 
1.000 
1.000 
0.998 
0.994 
0.992 
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U/Us 
1.011 
1.014 
1.018 
1.018 
1.013 
1.009 
1.011 
1.008 
1.009 
1.010 
1.007 
1.008 
1.003 
1.028 
1.065 
1.060 
1.028 
1.023 
1.020 
1.008 
1.006 
1.008 
1.009 
1.007 
1.004 
1.004 
1.002 
1.001 
1.001 
0.995 
0.993 


U/Us 
1.022 
1.012 
1.010 
1.011 
1.018 
1.010 
1.016 
1.012 
1.009 
1.009 
1.009 
1.005 
1.009 
1.105 
1.147 
1.184 
1.144 
1.092 
1.014 
1.008 
1.009 
1.006 
1.005 
1.008 
1.007 
1.007 
1.005 
1.001 
1.000 
0.996 
0.992 
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T% 
0.589 
0.681 
0.754 
1.229 
2.751 
4.241 
5.189 
3.855 
4.948 
4.298 
2.186 
3.431 
3.292 
6.959 
7.895 
9.544 

10.665 
9.366 
4.466 
2.118 
1.404 
2.895 
2.112 
2.824 
2.217 
2.517 
2.255 
1.471 
0.768 
0.584 
0.621 


U/Use 
1.018 
1.021 
1.037 
1.080 
1.068 
1.020 
0.990 
0.979 
0.938 
0.935 
0.995 
1.019 
1.006 
1.090 
1.142 
1.175 
1.175 
1.020 
1.004 
1.014 
1.006 
0.997 
1.013 
0.991 
1.016 
1.068 
1.079 
1.017 
0.997 
0.999 
0.993 
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U/Us 
1.019 
1.025 
1.042 
1.072 
1.058 
1.013 
0.934 
0.954 
1.025 
0.996 
0.966 
0.992 
1.051 
1.130 
1.130 
1.081 
1.065 
1.074 
1.012 
0.975 
0.982 
1.015 
1.046 
1.010 
0.978 
1.021 
1.067 
1.052 
0.999 
0.996 
0.993 


T% 
1.134 
2.154 
6.770 
6.771 
6.900 
5.519 
8.899 
10.237 
10.417 
12.672 
8.389 
7.359 
6.326 
7.105 
10.775 
9.992 
11.571 
8.133 
7.690 
7.178 
7.611 
8.822 
8.041 
11.003 
7.761 
9.185 
6.822 
7.477 
2.334 
1.031 
0.726 


U/Us 
1.018 
1.018 
1.022 
1.034 
1.038 
0.997 
0.969 
0.965 
1.098 
1.079 
1.129 
1.087 
1.118 
1.071 
1.035 
0.965 
0.969 
1.031 
1.062 
1.021 
1.080 
1.090 
1.001 
0.966 
0.953 
1.033 
1.067 
1.023 
1.005 
0.998 
0.993 


T% 
0.597 
0.966 
3.009 
3.560 
6.509 
8.447 

10.748 
14.355 
10.141 
10.349 
10.492 
11.327 
9.081 
9.011 
10.872 
9.726 
8.355 
10.691 
9.844 
10.789 
11.130 
11.761 
13.730 
10.879 
11.082 
8.672 
8.627 
6.134 
1.998 
1.379 
0.669 
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-2.0 


-2.0 
-2.0 


-2.0 
-2.0 
-2.0 


-2.0 
-2.0 


-2.0 
-2.0 
-2.0 


-2.0 
-2.0 


-2.0 
-2.0 
-2.0 


-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 
-2.0 


U/Us 
1.022 
1.021 
1.018 
1.022 
1.034 
1.030 
1.023 
1.035 
0.919 
0.914 
0.986 
0.981 
1.023 
0.968 
0.885 
0.873 
0.885 
0.868 
1.000 
1.018 
1.017 
0.966 
0.977 
1.013 
1.072 
1.075 
1.070 
1.007 
1.000 
0.998 
0.993 


T% 
0.622 
0.711 
1.369 
1.405 
4.712 
5.757 
8.176 

10.319 
13.444 
15.714 
12.408 
12.581 
11.645 
11.867 
9.502 
10.770 
14.755 
13.186 
12.909 
10.611 
13.019 
12.347 
10.748 
9.848 
8.903 
10.197 
9.761 
4.607 
1.777 
0.814 
0.689 


U/Us 
1.021 
1.022 
1.024 
1.022 
1.021 
1.020 
1.034 
1.075 
1.033 
0.960 
0.920 
0.923 
0.920 
0.929 
0.944 
0.956 
0.936 
0.927 
0.858 
0.880 
0.864 
0.973 
1.062 
1.076 
1.042 
1.033 
1.005 
1.008 
1.004 
0.998 
0.996 


T% 
0.651 
0.596 
0.727 
2.819 
2.376 
4.318 
6.015 
9.097 
9.717 

10.584 
9.807 
11.303 
11.386 
11.880 
12.118 
9.947 
8.700 
10.183 
10.433 
12.134 
9.239 
9.845 
10.686 
7.769 
7.660 
11.112 
8.037 
3.928 
1.175 
0.765 
0.738 
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-4.0 


-4.0 
-4.0 


-4.0 
-4.0 
-4.0 


-4.0 
-4.0 


-4.0 
-4.0 
-4.0 


-4.0 
-4.0 


-4.0 
-4.0 
-4.0 


-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 
-4.0 


U/Use 
1.023 
1.024 
1.028 
1.026 
1.026 
1.037 
1.039 
1.038 
1.035 
0.997 
0.961 
1.019 
1.043 
1.020 
1.029 
1.084 
1.028 
1.098 
0.983 
0.982 
0.958 
1.019 
1.025 
1.028 
1.016 
1.014 
1.013 
1.007 
1.009 
1.002 
0.997 


T% 
0.610 
0.538 
0.709 
1.691 
3.247 
3.051 
5.212 
7.747 
10.315 
12.606 
9.971 
12.892 
12.607 
10.869 
14.480 
12.333 
12.652 
12.667 
13.428 
12.231 
9.533 
10.298 
7.568 
8.688 
6.004 
3.220 
1.516 
1.076 
0.843 
0.681 
0.562 


U/Use 
1.034 
1.031 
1.035 
1.020 
1.032 
1.046 
1.108 
1.105 
1.050 
1.000 
1.080 
1.086 
1.151 
1.127 
1.152 
1.215 
1.222 
1.102 
1.046 
0.959 
1.001 
1.000 
1.053 
1.020 
1.065 
1.025 
1.008 
1.007 
1.007 
1.004 
1.000 


T% 
0.793 
0.751 
0.900 
2.407 
4.784 
7.950 
10.218 
11.943 
12.836 
12.657 
15.496 
12.165 
12.785 
12.514 
12.473 
12.944 
12.986 
14.956 
13.347 
12.983 
10.733 
12.729 
13.959 
17.636 
9.909 
4.880 
3.499 
1.225 
1.480 
0.708 
0.759 
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